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The  subject  contained  in  this  annual  report  covers  the  task 
"Statistical  Improvement  of  Ground-Based  Infrared  Profiles." 

During  the  past  year,  we  have  developed  a simple,  yet  sufficiently 
accurate  atmospheric  transmission  model  to  be  used  in  the 
15  jun  CO2  absorption  region  for  the  inversion  of  atmospheric 
temperature  profile.  The  result  of  this  study  has  been 
accepted  for  publication  in  June  issue  of  Optics  Letters.  A 
preprint  of  this  article  is  included  in  this  report.  Basically, 
we  Included  temperature  effect  of  the  transmission  function  in 
the  model  such  that  the  RIAPS  could  be  operated  with  seasonal 
as  well  as  geographical  variations.  If  a set.  of  transmission 
functions  is  known  either  analytically  or  empirically  for  a 
given  temperature  distribution,  the  first-order  approximation 
method  developed  here  may  be  used  to  calculate  the  transmission 
functions  for  a different  temperature  distribution.  Since 
coefficients  of  the  linear  function  are  predetermined,  this 
procedure  can  be  implemented  in  an  iteration  process  without 
significantly  Increasing  the  computation  time. 

The  new  transmission  model  is  yet  to  be  Incorporated  in  the 
RIAPS  software  program.  With  current  available  measurement 
data  (taken  in  the  summer  at  San  Diego  and  Pt.  Mugu) , it  is 
difficult  to  verify  effectiveness  of  the  new  transmission  model. 
However,  numerical  simulation  approach  may  be  employed  to 
compare  the  retrieved  profiles  using  existing  and  modified 
transmission  functions.  With  existing  line-by-line  calculation 
code,  one  can  further  look  into  combinations  of  various  frequency 
channels  with  distinct  transmission  characteristics.  The 
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discrepancy  between  calculated  and  measured  radiances  at 
certain  frequency  regions  is  a known  fact,  combining  RIAPS 
measurements  and  analytic  code,  one  might  be  able  to  shed  some 
light  on  this  problem. 
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Abstract 


The  method  of  first-order  Taylor  expansion  is -used  to  study  the  temperature 
effect  on  the  transmission  function  for  several  narrow  spectral  intervals  in  the 
CO2  absorption  region  with  direct  application  to  the  problem  of  thermal 
remote  sounding.  Numerical  analysis  indicates  that  the  computed  transmit- 
tances  are  in  good  agreement  with  those  obtained  by  line-by-line  calculation, 
especially  for  spectral  intervals  with  relatively  strong  temperature  depen- 
dence. — - 
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Infgnnation  on  the  transmission  characteristics  of  a real  atmosphere  in  narrow 
spectral  intervals  is  nece^ary  to  derive  atmospheric  temperature  and 
humidity  profiles  from  satellite  or  ground-based  radiation  measurements.  The 

16-ft  m CO2  absorption  band  has  been  chosen  for  atmospheric  thermal  remote  1 

sounding,  since  carbon  dioxide  is  uniformly  distributed  (assuming  330  ppm  in 

this  study)  through  the  entire  atmosphere.  In  a previous  study',  a set  of 

ground-based  irradiance  measurements  with  spectral  resolution  around  7 cm*' 

were  taken;  the  band  centers  are  located  at  669,  677.5,  688,  693,  700,  711.5, 

724,  and  730  cm*'.  These  measurements  were  related  to  the  atmospheric  tem- 
perature distribution  through  the  radiative  transfer  equation.  With  lineariza- 
tion of  the  Planck  function,  the  temperature  distribution  was  retrieved  through  i 

a linear  matrix  inversion  of  the  Fredholm  integral  equation  of  the  first  kind. 

An  iteration  process  was  included  to  correct  the  influence  of  the  water  vapor  j 

distribution  obtained  from  irradiance  measurements  in  the  18-^  m region. 

Since  this  inversion  problem  involves  an  iteration  process,  atmospheric  | 

transmission  functions  need  to  be  computationally  efficient  and  highly  accur-  | 

ate.  Originally  these  functions  were  empirically  determined,  assuming  tern-  \ 

perature  independence;  obviously,  this  treatment  of  the  transmission  functions  ! I 

I ^ 

is  not  adequate  under  conditions  of  substantial  seasonal  or  geographic  tem-  ; 

perature  variations.  One  way  to  improve  formulation  of  the  transmission  func- 
tion is  to  include  the  temperature  effect. 

There  are  numerous  methods  to  compute  the  atmospheric  transmittance 
and  the  most  accurate  one  is  line-by-line  calculation.*  In  this  procedure,  the 
transmittances  in  narrow  spectral  intervals  are  computed  as  sums  of  contribu- 
tions from  each  of  the  lines  in  the  absorption  band  and  are  integrated  over  the 
atmospheric  path.  This  approach  is  much  too  laborious  and  time  consuming  for  • 

use  in  real-time  processing  of  remotely  sensed  data. 

Alternative  approaches  include  the  Curtis-Godson  method  and  the  effec- 
tive mass  method*,  which  replace  the  inhomogeneous  path  with  an  equivalent 
homogeneous  path  by  defining  effective  values  of  pressure,  temperature,  and 
quantity  of  absorber.  This  approach  is  efficient  but  not  sufficiently  accurate. 

Recently  several  approximate  methods  have  been  proposed  for  fast  i 

calculation  of  the  atmospheric  transmittance  for  satellite  radiation  meausure- 
ments.  An  alternative  method  is  presented  here  using  transmission  functions  ' 

for  the  ground-based  measurements  in  the  numerical  analysis.  The  current  ap- 
proach, originally  suggested  by  Kondratyev  and  Timofeyev*,  expands  the 
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transimttance  at  a certain  pressure  level  in  Taylor  series  using  temperature  as 
the' variable.  This  method  is  especially  useful  when  the  transmission  functions 
for  a partictUar  temperature  distribution  are  accurately  known,  and  when  the 
approximate  method  is  used  to  compute  a correction  term  of  the  second-order 
of  smallness  for  a different  temperature  distribution.  Transmittance  between  x 
■■  0 and  X|  of  a nonhomogeneous  atmosphere  over  a spectral  interval  AP  may 
be  written  as 


‘»'a>  <«i>  - exp  u(x)  ^ Ki  (P,T,x)  dx]  d P / dP 


where  0(i^  is  the  instrument  response  function,  assumed  to  be  triangular  in 
this  analysis;  u is  the  optical  thickness  in  the  unit  of  atm -cm;  and  is  the  ab- 
sorption coefficient  of  the  ith  line.  In  the  lower  atmosphere,  the  Lorentz  line 
profile  for  the  absorption  coefficient  is  adequate: 


Kt  (p,T,x)  - 


a,Si(T) 


r((P-Pl)*  + ai*l 


The  line  intensity  SjCT)  for  a linear  molecule  may  be  written 


S,(T)  - S,(To)  (■“)*«P  [-  1.439  El  1 


where  T^  is  a reference  temperature,  Ei  is  the  lower  energy  level  of  the  absorp- 
tion line,  and  the  intensity  is  presented  in  the  unit  of  atm'>  cm'*  at  temperattire 
T.  The  induced  emission  and  the  vibrational  partition  function  have  been  ig- 
nored in  Eq.  (3).  The  line  half -width  ai  is  approximately  related  to  a reference 
value  ttiQ  at  pressure  and  temperature  T^  as: 


P 


The  temperature  dependence  of  the  half-width  is  neglected  in  Eq.  (4),  since  its 
influence  on  the  transmittance  is  much  smaller  than  that  due  to  the  line  inten- 
sity for  the  temperature  range  being  considered.  Substituting  Eqs.  (2)-(4)  into 
» Eq.  (1),  the  transmission  fiinction  is  expressed  explicitly  in  terms  of  the  tem- 
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The  term  Sr^p/  5T  is  the  variational  derivative  which  may  be  expressed 


Ki  [ -|^  + dv  / 0(1/)  d*/ . (7) 

From  Eqs.  (5)  through  (7),  the  transmission  function  for  a temperature  dis- 
tribution T(x)  may  be  estimated  when  the  transmittances  for  the  reference 
temperature  distribution  T^ix)  are  known.  For  a homogeneous  atmosphere,  the 
formulation  in  Eq.  (5)  is  reduced  to  the  familiar  form 

r^v  It^(T-T,).  (8) 

For  numerical  analysis,  th^  transmittance  T^p  [T(x)|  and  T^p  [To(x)],  and 
the  variational  derivative  dT^p/iiT  are  computed  using  the  line- by-line 
calculation  method.  The  approximate  transmittance  T^p  [T(x)]  is  then  evalu- 
ated from  Eq.  (5)  and  compared  with  the  value  obtained  from  the  line-by-line 
method.  The  line  parameters  were  generated  using  molecular  constants  con- 
tained in  the  report  of  McClatchey,  et  al*.  The  spectral  intervals  employed  in 
Ref.  1 were  chosen  to  investigate  the  accuracy  of  the  current  method.  The  in- 
terval at  669  cm*'  is  practically  opaque  for  surface  temperature  measurement, 
while  the  transmisison  function  at  724  cm'*  is  very  close  to  that  at  730  cm'*. 
These  two  channels  are  therefore  not  considered  in  the  current  analysis. 
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Since  t.he  temperature  dependence  of  the  transmiasion  function  is  determined 
chiefly  by  the  line  intensities  expressed  in  Eq.  (3).  lines  with  different  tern* 
perature  dependence  in  a spectral' interval  could  lead  to  difflculties  in  deter* 
mining  the  dependence  of  the  transmission  function  on  temperature.  The 

nature  of  this  problem  may  be  seen  from  Eq.  (7),  where  the  two  terms-nrand 
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t: — are  opposite  in  sign.  Thus,  for  a spectral  interval  where  most  of  the 

T* 

lower  energy  levels  of  strong  absorption  lines  are  greater  than  a certain  value 
(•«1.39T),  the  transmittance  increases  as  the  temperature  increases  within  a 
certain  temperature  range.  For  the  range  220K  to  340K,  the  transmittance 
decreases  as  the  temperature  increases  for  the  677.5  cm*>  spectral  interval, 
while  other  channels  show  the  opposite  trend.  The  temperature  dependence  at 
688  cm'>  is  very  weak;  the  first -order  approximation  method  failed  to  improve 
the  transmittance,  since  a much  more  accurate  numerical  integration  scheme 
is  required.  Figure  1 is  a plot  of  the  transmittances  as  a function  of  the  optical 
thickness  for  a homogeneous  atmosphere  at  300K  and  1 atm.  These  results  were 
obtained  by  line-by-line  calculations  and  first-order  approximation  using  two 
reference  temperatures,  260K  and'340K.  Usually,  the  difference  between  the 
reference  temperatxire  and  the  actual  temperature  is  less  than  40K  (chosen  in 
this  example).  Apparently  use  of  a higher  reference  temperature  gives  better 
agreement  between  the  approximate  transmittance  and  the  ideal  transmit- 
tance. Furthermore,  agreement  between  actual  and  approximate  transmit- 
tances is  always  better  in  an  optically  thin  atmosphere.  The  pereentage  of  error 
of  the  approximate  method  is  plotted  in  Figure  2.  For  those  portions  of  the 
curves  with  large  errors,  the  actual  transmittance  is  less  than  0.001;  thus,  the 
comparison  is  physically  less  significant. 

This  approximate  method  is  also  evaluated  for  a nonhomogeneous  at- 
mosphere. First,  the  transmission  fimctions  for  the  1962  standard  atmosphere 
and  the  tropical  model  atmosphere^  are  computed  using  line-by-line  calcula- 
tions, and  then  the  approximate  transmission  fimctions  for  the  1962  standard 
atmosphere  are  estimated  using  the  tropical  atmosphere  as  the  reference  tern- 
perture  distribution.  These  results  are  listed  in  Table  I where  transmittance 
values  less  than  0.001  are  entered  as  zero.  In  most  cases,  the  accuracy  of  the 
approximate  method  is  within  one  percent  of  the  actual  value,  and  the  worst 
discrepancy  is  within  10  percent,  which  occurs  at  677.5  cm'*.  The  approximate 
transmittances  at  688  cm'*  were  not  computed,  since  a more  accurate  numeri- 
cal integration  scheme  is  required  for  this  channel. 


In  summary,  the  temperature  dependence  of  the  transmission  function  for 
narrow  spectral  intervals  in  the  15-m  m region  has  been  investigated  using  the 
method  of  Taylor’s  first'order  approximation.  Except  for  688  cm*>,  the  com* 
puted  transmittance  using  the  approximate  method  is  in  good  agreement  with 
that  using  the  line-by-line  calculation  method  when  the  reference  temperature 
distribution  is  properly  aosen.  This  method  is  highly  accurate  and  computa- 
tionally efficient  for  application  to  thermal  remote  sensing  problems.  It  is 
especially  suitable  when  transmission  functions  are  accurately  known  for  a 
specific  model  atmosphere  and  corrections  are  computed  using  the  current 
method.  Although  the  numerical  analysis  presented  here  is  for  a ground-based 
thermal  sounding,  it  should  apply  equally  to  the  problem  of  downward-looking 
satellite  radiation  measurements. 

This  work  was  supported  by  the  Office  of  Naval  Research  under  contract 
number  N00014-75-C-0940. 
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Ft^re  1.  Comparison  of  computed  COy  transmittanees  using  line-by-line  caleula-  j 

tion  and  first-order  approximations  spectral  resolution  is  7 cm'K 


Figure  2.  Percentage  of  error  of  the  computed  transmittanees  using  the  approxi- 
mate method  for  the  same  spectral  intervaU  as  in  Figure  1. 


Thble  L COj  transmittanees  between  the  ground  and  selected  altitudes  for  several 
spectral  intervals  in  the  15-p  m region. 
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Table  L CO2  transmittances  between  the  ground  and  selected  altitudes  for  several 
spectral  intervals  in  the  IS-p.  m region. 


Frequency 

(cm*^ 

Approach* 

Temperature 

Profile** 

Altitude  (km) 

0.1 

0.3 

0.5 

0.8 

1.5 

3.0 

6.0 

877.5 

L 

S 

0.108 

0.0047 

0 

0 

0 

0 

0 

F 

S 

0.112 

0.0051 

0 

0 

0 

0 

0 

L 

T 

0.119 

a0059 

0 

0 

0 

0 

0 

' 688 

L 

S 

0.161 

0.0215 

0.0042 

0 

0 

0 

0 

L 

T 

0.164 

0.0225 

0.0045 

0 

0 

0 

0 

693 

L 

S 

0.253 

0.0480 

0.0117 

0.0019 

0 

0 

0 

F 

S 

0.258 

a0497 

0.0123 

0.0020 

0 

0 

0 

L 

T 

0.247 

0.0460 

0.0111 

0.0017 

0 

0 

0 

700 

L 

S 

0.396 

0.133 

0.0536 

0.0163 

0.0017 

0 

0 

F 

S 

0.398 

0.135 

0.0544 

0.0166 

0.0017 

0 

0 

L 

T 

0.377 

0.119 

0.0456 

0.0131 

0.0012 

0 

0 

711.5 

L 

S 

0.713 

0.475 

0.357 

0.251 

0.135 

0.0584 

0.0277 

F 

S 

0.712 

0.477 

0.359 

0.254 

0.136 

0.0587 

0.0276 

L 

T 

0.690 

0.450 

0.330 

0.226 

0.115 

0.0460 

0.0200 

730 

L 

S 

a866 

0.718 

0.633 

0.553 

0.449 

0.342 

0.263 

F 

s 

0.867 

0.720 

0.635 

0.557 

0:452 

0.347 

0.267 

L 

T 

0.860 

0.709 

0.622 

0.542 

0.436 

0.327 

0.245 

* L - line-by-line  calculation  method;  F - first-order  approximation 
^ S - 1962  standard  atmosphere;  T - tropical  model  atmosphere 


